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ABSTRACT 

Some galaxy clusters show diffuse radio emission in the form of giant halos (GHs) on Mpc 
scales or minihalos (MHs) on smaller scales. Comparing Very Large Array and XMM-Newton ra- 
dial profiles of several such clusters, we find a universal linear correlation between radio and X- 
ray surface brightness, valid in both types of halos. It implies a halo central emissivity vj^ — 
10-3i-4±0-2(^/XO-2 cm-3)2(T/Tn)°-2±0-5 erg s'^ cm.-^, where T and Tq are the local and central tem- 
peratures, and n is the electron number density. We argue that the tight correlation and the scaling 
of combined with morphological and spectral evidence, indicate that both GHs and MHs arise 
from secondary electrons and positrons, produced in cosmic-ray ion (CRI) collisions with a strongly 
magnetized, B > 3 intracluster gas. When the magnetic energy density drops below that of the 
microwave background, the radio emission weakens considerably, producing halos with a clumpy mor- 
phology {e.g., RXC J2003. 5-2323 and A2255) or a distinct radial break. We thus measure a magnetic 
field B — 3 /iG at a radius r ~ 110 kpc in A2029 and r ~ 50 kpc in Perseus. The spectrum of 
secondaries, produced from hadronic collisions of ~ 20 GeV CRIs, reflects the energy dependence of 
the collision cross section. We use the observed spectra of halos, in particular where they steepen 
with increasing radius or frequency, to (i) measure B ~ 10(i//700 MHz) /iG, with i/ the spectral break 
frequency; (ii) identify a correlation between the average spectrum and the central magnetic field; and 
(in) infer a CRI spectral index s < —2.7 and energy fraction ~ ]^q-3.6±o.2 particle energies above 
10 GeV. Our results favor a model where CRIs diffuse away from their sources (which are probably 
supernovae, according to a preliminary correlation with star formation), whereas the magnetic fields 
are generated by mergers in GHs and by core sloshing in MHs. 

Subject headings: galaxies: clusters: general — intergalactic medium — X-rays: galaxies: clusters — 
radio continuum: general — magnetic fields 



1. INTRODUCTION 

Giant halos (GHs) appear as diffuse radio emission 
on ~ M pc scales in merging galaxy clusters (for a re- 
view, see lFeretti fc Giovanninil 120081) . GHs were identi- 
fied in about a quarter of all clusters with X-ray lumi- 
nosities Lx > 5 x 10^ ** erg s~^ at redshifts 0.2 < z < 0.4 
(|Brunetti et al.|[2007D . Lx and the specific radio power 
Pu of GH clusters are tightly correlated. The GH distri- 
bution in clusters is bimodal, with most clusters showing 
no associated GH at a sensitivity threshold ^ 10 times 
bette r than the signal exp ected from the Lx-Pu correla- 
tion (iBrunetti et al.ll2007[ ). 

GHs arise from synchrotron radiation, emitted by 
cosmic-ray electrons and positrons (CREs) injected lo- 
cally and continuously into the magnetized plasma. 
Such CREs can lose a considerable fraction of their 
energy by inverse-Compton scattering off the cosmic 
microwave background (CMB). Recentlv. iKushnir et al.l 
(|2009l henceforth K09) made the important observation 
that the GH properties mentioned above are reproduced 
if the CREs arc produced in p-p cosmic ray proton (CRP) 
collisions, and the magnetic field B is sufficiently strong 
to saturate Pu{B) by rendering Compton losses negligi- 
ble. 

To qualitatively see this, assume that the CRP num- 
ber density Up is (narrowly distributed about) a univer- 

^ Einstein fellow 



sal fraction of the local (non cosmic-ray) electron number 
density n, and that the magnetic energy density greatly 
exceeds the energy density Ucmb of the CMB, B ^ 
Bcmb = (STTUcrnf,)^/^. Here, the ratio between the emis- 
sivities of radio synchrotron jy oc UpuB^ / {B^ + B^^^f^) cx 
(from secondary, density ne oc npU CREs) and X- 
ray bremsstrahlung jx oc (from the thermal plasma), 
does not depend on n or on B. Thus, strongly magne- 
tized clusters with B > Bcmb would show a strong radio- 
X-ray correlation, whereas clusters with B < Bcmb/i 
would be too radio faint to show a GH. Other models, 
not ably turbulent re acceleration of electrons (for a review 
see iPetrosian fc Bvk ov 2008), do not naturally produce 
the correlation and bimodality that are observec0. 

Radio minihalos (MHs) are found in cool core clus- 
ters (CCs). They extend roughly over the cooling region 
(jGitti et al.ll2002t ). encompassing up to a few percent of 
the typical GH volume. Detecting MHs is more challeng- 
ing than GHs, due to their smaller size and proximity to 
an active galactic nucleus ( AGN) , so only few MHs have 
been well studied and less is known about their correla- 
tion with X-rays. A morpholo gical association between 
MH e dges and cold fronts (CFs: lMazzotta &: Giacintuccil 
120081) suggests a link between MHs and sloshing ac- 
tivity in the core. Such CFs are observed in about 

In the sense that the small dispersi on above PyjLx) is not re - 
produced, and many assumptions (e.o..'Brunetti fc Lazarianl l2007l ) 
are needed. For a different opinion, see Brunctti et al. ( 20091) . 
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half the CCs, and are probab ly present in many more 
(|Markevitch fc Vikhlininll2007[ ). They are thought to be 
tangential discontinuities that isolate regi ons magnetized 
by bulk shear flow at smaller radii (Kes het et al.l[20Tol ) 
(regions we refer to as below, or inside the CF). 

There are many similarities between GHs and MHs. 
Both types of halos are usually characterized by a reg- 
ular morphology, low surface brightness, little or no po- 
larization, and spectral indices = (ilog(Pi,)/(ilogi' ~ 
— (1.0 — 1.5), with V the frequency. However, there are 
telling exceptions to these characteristics. A few GHs 
have a clumpy or filamentary morpholog y, such as in 
RXCJ 2003. 5-2323 (iGiadntucci et al.l l2"00l. A2255, and 
A2319 (|Murgia et al.ll2009l henceforth M09). Strong po- 
larization was so far dete cted in one GH (at a 20%- 
40% level, in A2255; see iGovoni et all (|2005t ): weaker 
polarization, 2 — 7 % on average, w as found in MACS 
J0717.5 -h3745; see iBonafede etall (i200a) ), and in one 
MH (at a 10% - 20% level, in A2390; iBacchi et al] 
|2003() . Spectral steepening with increasing r, increas- 
ing V, or decreasing T, has been reported in several ra- 
dio halos CFeretti fc Gi ovanninill200a iFerrari et alll2008l : 
iGiovannini et al. 2009, and references therein); typically, 
a steepens from —1.0 to —(1.3-1.5) in uncontaminated 
regions (see tj5.2p . A handful of GHs show a steep, 
a < —1.5 spectrum, s uch as in A521 {a = —1.86 ± 0.08; 
Dallacasa et al.' '2009') and in A697 (a = -(1.7-1.8); 
Macario et al. 2010). 

We analyze a sample of GHs and MHs using radio and 
X-ray data from the literature. We find a universal cor- 
relation between the radio and X-ray surface brightness, 
which holds for both types of halos. This correlation, 
combined for example with the observed dependence of 
the radio bright volume fraction upon cluster parame- 
ters fCassano et al. 2007), gives rise to the luminosity 
correlation known in GHs, and a similar correlation that 
we derive for MHs. We determine the radio emissivity 
ji, and its scaling with electron density and temperature 
T, using the surface brightness correlation and a model 
(XSPEC/MEKAL) for the X-ray emission. Combined 
with other properties of GHs and MHs, this singles out 
secondary CRE models with strong magnetic fields for 
both types of halos. This generalizes the GH model of 
K09, and, considering the different environments of GHs 
and MHs, substantially strengthens it. We propose that 
the cosmic-ray ions (CRIs) are produced in supernovae, 
based on the halo spectra, the ji, scaling, and a prelimi- 
nary correlation between star formation and the radio to 
X-ray brightness ratio rj. We show how the spectral and 
morphological properties of halos can be used to disen- 
tangle the distributions of CREs and magnetic fields. In 
particular, the magnetic field is gauged both by a radial 
break in rj reflecting the onset of Compton losses, and by 
a spectral break in radio emission induced by diffractive 
p-p scattering. 

The paper is arranged as follows. In fJ5] we show that 
GHs and MHs follow the same correlation between spe- 
cific radio power and coincident X-ray luminosity, sug- 
gesting that the different types of halos arise from similar 
processes. A tight, universal correlation between the sur- 
face brightness in radio and in X-rays is presented in 
It is used to derive the scaling of radio emissivity with 
the local electron number density n and temperature T, 
and to explain previous results for the average halo emis- 



sivity and the luminosity correlation. We then show in 
21 that the inferred scaling favors a model in which the 
radio emission in both GHs and MHs originates from 
secondary CREs produced by collisions of CRIs with the 
intracluster gas, synchrotron radiating most of their en- 
ergy in strong magnetic fields. We study the morpholog- 
ical and spectral properties of the halos in SjSl presenting 
new methods for measuring the distributions of CRIs and 
magnetic fields. Finally, we discuss the origin of the CRIs 
and magnetic fields in ^ and summarize our results. We 
assume a Hubble constant H — 70 km s~^ Mpc~^. Error 
bars are la confidence intervals. 

2. GHS AND MHS: SIMILAR CORRELATION BETWEEN 
RADIO AND COINCIDENT X-RAY EMISSION 

The 19 X-ray bright GH clusters in the lBrunetti et all 
(j2007| ) sample follow the correlation 

where the subscript 45 denotes units of 10^^ ergs~^. 
Here, X[ei,e2] designates integration over photon en- 
ergies between ei and €2 measured in keV (we use 
X[0. 1,2.4] unless otherwise stated), and subscript 1.4 
means that the radio signal is evaluated at a frequency 
v = \ A GHz. The normalization uncertainty in Eq. ([T]) 
is dominated by an intrinsic scatter among halos, needed 
to obtain an acceptable fit (chosen as /N — 1 with N 
being the number of degrees of freedom). Accounting for 
dispersion in P^, leads to the correlation in the form of 
E g. dH), in agree r uent w ith K09, somewhat flatter than 
m (Cassano et all (|2006[ ). 

A sample of 6 we ll-studied MHs had been reported 
(jCassano et al."2OO80 as being "barel y " con sistent with 
the GH correlation of ICassano et al.l (|2006[ ). However, 
we find that this sample agrees well with Eq. ([1]) with 
X^/N = 1.4. Nevertheless, within a partly overlapping 
sample of 6 MHs studied recently (|Govoni et al.l 120091 
M09), the MHs appear significantly underluminous with 
respect to the prediction of Eq. (|T]) after a careful removal 
(M09) of the contamination from AGNs. 

These results do not imply that MHs are intrinsically 
fainter than GHs. Rather, MHs have comparable radio 
power, if one corrects for their smaller size. A natural 
generalization of Eq. ([T]) applicable for both GHs and 
MHs would replace Lx with the luminosity Lx of the 
radio bright region. Following M09, we define the radio 
bright region based on the surface brightness condition 
I^{r)/Iy{r = 0) > ~ 5%, with r being the radial 
distance from the cluster center; equivalently, the radius 
of this region is — Zr^, with Tg being the e-fold radius. 

A quick way to estimate Lx is to correct Lx using a 
model for the X-ray brightn ess profile of the cluster. Us- 
ing an isothermal /3-model (|Cavaliere fc Fusco-Femianol 
1 1976( 1. in which n(r) = no(l -I- r^/r^)"^''/^ for model pa- 
rameters fi, Vc and no, we find that 

r J + r2i 

Lx ^ Lx . (2) 

We use individual cluster /3-model parameters from the 
literature; see Table |T] for details. The radius Rx of the 
X-ray region considered must be finite in order to ensure 
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convergence. Its choice is somewhat arbitrary; we use 
Rx ~ i2rc so Ri, < Rx for ah halos in the M09 sample. 

Figure [T] shows various GHs and MHs in the result- 
ing Lx — vPv plane. Replacing Lx by Lx results in 
slightly better agreement of the M09 GHs with Eq. ([1]), 
yielding /N = 1.3 instead of 1.6 (before propagating 
/3 model uncertainties). However, the MH agreement be- 
comes much better, /N = 2.4 instead of 3.7. Note 
that the fit is not expected to be as good for MHs as it is 
for GHs, because a /3-model is less appropriate for a CC. 
These results suggest that the relation between radio and 
X-ray emission is similar in GHs and in MHs. For more 
accurate, model-independent results, we next examine a 
more local measure of the emission, and consider the ra- 
dio and X-ray morphologies. 
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Fig. 1. — Radio power vP^ at 1.4 GHz of GHs (red trian- 
gles) and MHs (black circles) from M09, plotted against the 0.1- 
2.4 keV X-ray luminosity Lx of the radio bright region (see 
ij2] and table [T] for details). The total X-ray luminosity Lx is 
used in cases where Lx could no t be determined, for GHs (blue 
squares, from 'Cassano et al.!l2fl06!) and MHs ( green diamonds, from 
ICassano et a l. 2008). The solid line and shaded region show the 
best fit and la interval of the GH correlation in Eq. JlJ. 



3. CORRELATION BETWEEN RADIO AND X-RAY 
SURFACE BRIGHTNESS: UNIVERSAL EMISSIVITY 

A more useful manifestation of the radio- X-ray correla- 
tion is the morphological similarity between radio and X- 
ray emission. A linear correlation between the radio sur- 
fac e brightness / t, and th e X-ray brightness Fx was found 
bv lGovoni et all ()2001aD in two GH clusters (A2744 and 
A2255, inspected individually), while a sublinear power 
law was found in two other GHs (A2319 and Coma). 

We investigate the connection between the radio and 
the X-ray surface brightness in GH and MH clusters by 
examining the radial brightness profiles published in the 
literature. We consider not only the radio to X-ray re- 
lation within each cluster, but also compare the ratio 
between vl^, and in different clusters. Thus, Figure [2] 
shows the radial profile of the dimensionless ratio 

_ ^/.(1.4 GHz) 

= — (3) 

^X[0.1-2.4] 

between radio and X-ray surface brightness, for all six 
GHs and four MHs with published radial profiles from 
both Very Large Array (VLA; M09) and XMM-Newton 



(jSnowden et al.|[2008D observations. Distances are nor- 
malized to rsoo, the radius enclosing 500 times the criti- 
cal de nsity of the Un iverse (calculated using the best fit 
of . Zhang et al.l[2008h . 
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Fig. 2. — Radial profiles of the radio-to-X-ray ratio rj (defined 
in Eq. l(3}) in the GHs (open symbols) and MHs (filled symbols) 
with overlapping pro files from VLA (M09) and from XMM-Newton 
dSnowden et al.ll2008l . k-corrected). The MH data is shown only in 
regions where the AGN contamination is < 10%; this is supple- 
mented by the exponential model fits (dashed curves) given by 
M09, plotted in the range where radio detection is > Scr and the 
MH contribution is > 30%. Error bars are the la confidence inter- 
vals of the radio normalization /q (M09), and the solid lines serve 
to guide the eye. The la best fit of Eq. ^ is shown as a thin 
dashed line in a shaded reg ion. Model fits for A2029 (dotted and 
dot-dashed curves; see ilS.lll are also shown. 

All halo clusters except the MH in A2029 converge on 
a similar value of rj at small radii. The innermost GH 
data is best fit by 

= il{r = 0) = 10-4-4±o.2 ^ 

where the uncertainty is again larger than the measure- 
ment errors and mostly attributed to intrinsic scatter 
between halos. The inner MH data also lie within the 
la scatter, except A2029 which is significantly fainter 
in radio (but is observed only where ri(r) rapidly de- 
clines). Therefore, Eq. ^ appears to hold at small radii 
for both GHs and MHs. A possible explanation for the 
relatively high radio brightness of Perseus and A665, and 
the low brightness of A2029, is discussed in ij6.2l As in 
we used only GHs, which are better understood than 
MHs, to derive the correlation Eq. (jH). If, instead, we 
use the innermost data of both GHs and MHs, we obtain 
a broader best fit, ryo — 10^4.4±o.4^ ^j-^g bright MH 

in Perseus and the faint MH in A2029. 

Although there is some scatter among different clus- 
ters. Fig. [5] shows that rj is fairly uniform within each 
GH. An exception is A2163, in which rj{r) monotonically 
increases out to r ~ O.Brsoo — 1200 kpc. However, this 
is the only cluster show n in Fig. [2] that harbors a radio 
relic (iFeretti et al.ll200l[ ). which dominates the emission 
near the maximum (M09, figure 3). This suggests that 
the peculiar 77 profile in A2163 is associated with the 
radio relic. A peculiar 7y(r) pattern is seen in A2319, 
with a minimum at r ~ O.lrsoo ~ 150 kpc. Interestingly, 
this cluster shows both GH and MH characteristics, har- 
boring both a CF at r ~ 100-200 kpc and a subcluster 
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TABLE 1 

Parameters of Halo Clusters in the Sample 
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Columns: (1) cluster name (IRXS J131423 .6-251521 abbreviated RXJ1314; RXJ1347.5-1145 abbrev. RXJ1347); (2) halo type (GH or MH); (3) redshift z; (4) X-ray luminosity 
between 0.1 and 2.4 keV, in units of 1044 erg s~l; (5) temperature in keV; (6) e-fold radius of radio brightness in kpc (M09); (7) average 1.4 GHz specific emissivity in units of 
10~ erg s^i cm" -'■^ Hz~l (M09); (8-9) exponent /3 and core radius rc in kpc for an isothermal /3-model of the cluster; (10) X-ray luminosity of the radio bright region (see fj2)l between 
0.1 and 2.4 keV, in units of 1044 erg s~l; (11) minimal central magnetic field Bo^rnin iu units of /xG, assuming the /3-model and that B{r = minfr;,, 3re}) = i?cm6 (corresponding 
to fic = 0.5 in Eg. I12I I. where the radial break radius ri, = 55 kpc in Perseus, 135 kpc in A2029, 60 kpc in A2319, and ri, > 3re in the other halos; (12) average spectral index a'^^ 
between frequencies (in GHz) ui and 1^21 or olu around frequency v\ (13) Metallicity measured at r = O.lrsoo in solar units l lSnowden et al.ll2008l ). 

References: References arc given in parentheses. For conflicting references, we adopt the tighter estimate if applicable, and otherwise use the most recent result. AOS: 
[Allen ct al. (2003); A06: Augusto ct al. (2006); A09a: Andcrsson ct al. (2009); A09b: Anderson ct al. (2009); BOO: Bohringcr ct al. (20 00); B03: Bacchi ct al. (2003); B06: 
[Bonamcntc ct al. (2006); B07: Baldi ct al. (2007); C06: Cassano ct al. (2006); C07: ^Chen et al. (2007); D93: David ct al. (1993); F04: Fukaziwa et al.' (2004); F04b: Fcrctti ct J] 
(2004a) F07: Fcrctti ct al. (1997); GOl: Govoni ct al. (2001b); G03: Giovannini ct al. (2003); G04: Gitti ct al. (2004); H80: Harris ct al. (1980); KOI: KcmDner & Sarazin (200ll); 
MOO: .Matsumoto et al. (.2000 '): MOl: Matsumoto et a,!., (2001 1; M07: .Morandi ct al.. (2007 1; R81: Roland ct al (1981 ) R02: Reiprich k. Bohringcr 12002,'): S8 3: ,Slce fc Siegman (.19831 1: 
SOS: ISantos et al.l 1I2OO8II ; S09: [Sanderson et al.( (1200911 ; V02: IValtchanov et al.l (120031 : V09: Ivan Weeren et al.l (120091 1: WB03: IWorrall fc Birkinshawl (lioOSl I: WOO: IWhitel (I2000IV 
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at r ~ 500 kpc (|Govoni et al.ll2004[ ). In we argue 
that the clumpy radio morphology of this cluster sug- 
gests weak magnetization, enhanced at both small and 
large radii in an intermediate phase between a GH and a 
MH. The declining ri{r) profiles of the MHs are discussed 

in gm 

In order to derive the radio emissivity from the bright- 
ness correlation Eq. (|4]), we must model the X-ray 
emission. The X-ray emissivity calculated usin g the 
MEKAL m odel (M ewe et all [19851 [1986'; 'Kaastra 1992"; 
iLiedahl et al. 1995) in XSPEC v.12.5 (Arnaud 1996) is 
nearly independent of temperature and metallicity. It is 
given by 

^2 ^0.04 

JX[0.1-2A] ^ 8.6 X 10^28 -2 0.3 ^-1 ^^-3 ^ (5) 

where n^2 is the electron number density n in units of 
10-2 cm-3, Tio = (fcsT/lO keV) with ks being Bohz- 
mann's constant, and Zq 3 is the metallicity Z in units 
of O.SZq. 

Fairly uniform values of r/ are found within each cluster 
(excluding the outer regions of A2163), and similar val- 
ues are found among different halos, in regions spanning 
two orders of magnitude in density. This implies that 
the radio emissivity ji, , like jx , scales as the gas density 
squared. Indeed, M09 recently found (in their section 
5.1) that assuming j^, oc reproduces the radial profiles 
of some GHs and MHs. We conclude that the linear sur- 
face brightness relation Eq. ^ reflects a similar, local 
relation between the radio and X-ray emissivities, 

i^ji^ - Vojx ■ (6) 

We examine the temperature dependence of the radio 
emissivity within a given cluster by fitting 77 in uncontam- 
inated GHs as a power law, such that ji, cx n'^{T/To)^ , 
with To being the central temperature and k being a free 
parameter. This yields k = 0.2 ± 0.5, consistent with 
weak or no temperature dependence within the cluster. 
Combining this with the above results Eqs. dl])-®, we 
conclude that the synchrotron emissivity in halos is well 
fit by 

~ 10-31-4±0.2„^^(j./j.^)0.2±0.5 g-l ^^-3 ^ 

at least near the centers of the halos. 

The scatter in Eqs. @ and ([7]) indicates that the radio 
emissivity is not a function of the local density and tem- 
perature alone. The dispersion in probably reflects 
dependence upon additional cluster properties, such as 
the cluster's mass and star formation history. It has 
standard deviation cr(logji,) ~ log(1.6) among the GHs 
in our sample, or a{\ogj^) ~ log(2.5) if we include the 
MHs. This relatively small dispersion explains the simi- 
lar volume-averaged emissivity (j^) of different GHs, and 
the high level of variation in (j^) among MHs, where (ji,) 
is sometimes t wo orders of ma g nitud e larger than in GHs, 
as reported bv lCassano et al.l ()2008l ) and M09; see Table 
[TJ Indeed, the GHs in their sample have comparable sizes 
and central densities no — 10-^ cm"^, whereas the MHs 
have variable scale and uq > 10~^ cm"'^ in some cases. 

We have derived linear relations between the radio and 
X-ray emissivity (Eq. ^ and surface brightness (Eqs. [3| 
and[4]) of halos, but a superlinear, P^, cx L]. '^ relation be- 
tween the integrated luminosity in the two bands (Eq.[T]). 



This indicates that, in addition to the ji, cx n'^{T/To)'^ 
dependence, the total radio power of a halo must further 
increase, on average, with the mass and temperature of 
the cluster. Such a behavior could result in part from a 
(direct or indirect) mass or temperature dependence of 
ji,. However, the dispersion in is small, and we find 
no evidence for such a correlation. The effect is proba- 
bly dominated by the different scalings of the radio and 
X-ray bright volumes with the cluster parameters, as we 
qualitatively show next. 

The length scale i?^ of GHs was found to depend su- 
perlinearly on t he virial radius Ryir , Rv ^ ^tS-j with 
Xr = 2.6 ± 0.5 (jCassano et al.ll2007[ ). In order to esti- 
mate the effect of this scaling on the P^-Lx correlation, 
we crudely approximate the luminosity within a sphere 
of radius i? as a power-law, L{< R) cx R'^^ . In a /3- 
model, this is a poor approximation, with power-law in- 
dices 0.5 < Al < 2 in the parameter range 0.5 < (3 < 0.7, 
1 < r/rc < 4, relevant to GHs. For simplicity, we assume 
that Lx reflects the emission within a sphere of radius 
cx Rvir- For constant 77, we may then write 

Pu ( Rv \ r,(XR-l)\L j^0.3,(\r-1)\l fa\ 

— CX — CX Rl„ ' (xLx' , (8) 

where we used th e observed scalii igs Ryjr ^ T°-^ 
(jZhang et all [200g) and Lx ~ IjMarkevitchI [1991 
in the last proportionality (recall that X represents the 
energy range 0.1-2.4 keV). Equation ^ yields Pi, cx 
2^(1.2 2.0)^ where we accounted for the Al range but have 
not propagated the scaling errors. While this demon- 
strates that the R,^-Rmr scaling may sufhce to recon- 
cile the nonlinear luminosity relation Eq. ([T]) with the 
linear emissivity relation Eq. (|6]), much better modeling 
is required in order to identify all factors governing the 
luminosity correlation. For example, we neglected halo 
asymmetry and possible variations in at large radii (as 
in A2163). 

4. UNIVERSAL RADIO MECHANISM 

In the preceding sections we examined the luminosity 
and surface brightness properties of GHs and MHs, in ra- 
dio and in X-rays, and derived the radio emissivity, with- 
out making any assumptions regarding the particles and 
magnetic fields responsible for the radio emission. Here 
we explore the implications of the radio-X-ray correlation 
and the radio scalings derived above. In i i4.1[ we discuss 
various models for radio halos, and show that only mod- 
els with strong magnetic fields and secondary CREs are 
consistent with the observations. In M.2I we briefly dis- 
cuss the effect of time dependent magnetic fields. 

4.1. Secondary CREs, strong magnetic fields 

The spectral slope near the center of most radio ha- 
los is a ~ — 1 (see ^5.2p . A synchrotron spectrum of 
this type is emitted by rapidly cooling CREs, injected 
with approximately constant energy per decade in par- 
ticle energy E,, . The synchrotron emissivity correspond- 
ing to a logarithmic CRE energy density injection rate 
Q = due/{dt din Ee) = constant is given by 

. ^ du^ _ Q/2 
"^"-"dtdi^- l + iB,„.b/Br ' 
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where the denominator accounts for Compton losses, 
Bcmb = iS-KUcmbY^^ — 3.2(1 + z)^ is the amplitude 
of the magnetic field which has the same energy density 
as the CMB, and Uj is the energy density of component j 
(7 for radio photons, e for CREs, etc.). In the B Bcmb 
limit, the synchrotron emissivity Eq. ([7]) extracted from 
the observations yields a direct estimate of the CRE in- 
jection rate. 



Q 



Q 



1 + {Bcmb/Bf 



^Q-3i.i±o.2^2^_^ erg s"^ cni"^ 



(10) 



As CREs do not have time to diffuse before they lose 
their energy, this result applies locally. It holds in both 
GHs and MHs. 

The remarkably small scatter found in the radio-X-ray 
correlations among different GHs has been interpreted as 
implying a robust radio mechanism that keeps the syn- 
chrotron emissivity narrowly distributed around a uni- 
versal function of the plasma density and perhaps also 
temperature. Equation ([T0|) and the narrow scatter in 
its normalization suggest that the universal quantity in 
radio halos — both GHs and MHs — is Qjr? . 

Our results strongly emphasize the robustness of the 
radio emission mechanism. First, we find a tight radio- 
X-ray correlation not only in GHs but also in MHs. As 
the latter have physical properties that are considerably 
different from GHs (smaller size, higher ambient density, 
nearby AGN, CF association), the radio emission mech- 
anism must be sufficiently robust to reproduce the same 
levels of Qjr? over a wide range of physical conditions. 
Second, we confirm that a tight radio- X-ray correlation 
exists not only in the total cluster luminosity but also in 
the local surface brightness, and that the ratio t] between 
radio and X-ray brightness is fairly uniform within each 
cluster. This is particularly striking in MHs, where the 
coincident density and temperature profiles are steep. 

K09 have pointed out that in the strong magnetiza- 
tion regime, defined as i? > Bcmb, radio emission from 
GHs is independent of the precise value of B, and so the 
tight correlation only constrains the CRE injection Q. In 
contrast, in the weak magnetization regime B < Bcmb, 
the product QB^ should be universal, requiring a physi- 
cal mechanism that carefully balances CRE injection and 
magnetic fields both with each other and with the am- 
bient gas. Furthermore, in a strongly magnetized halo 
model, the transition from a constant to a rapidly de- 
clining (cx B^) behavior as B drops below Bcmb naturally 
explains the GH bimodality observed (K09). 

The similar values of ij we find in GHs and in MHs 
strengthen this argument considerably, because it is dif- 
ficult to come up with a double feedback mechanism 
(CRE-magnetic fields-ambient plasma) that operates 
identically in the different environments of GHs and 
MHs, without fine tuning. The bimodality argument 
does not apply to MHs, however, as their distribution 
has not been shown to be bimodal, and may in fact be 
continuous if high magnetization is ubiquitous in CC cen- 
ters (see discussion in pS.ip . 

Two types of models have been proposed for CRE in- 
jection: (i) seco ndary production by hadronic collisions 
involving CRIs (jPenniso nl [19801) : and (li) in-situ tur- 
bulent acceleration or reacceleration of primary CREs 
(|Enfilin et al.l |1999|) . These models typically assume 



fixed ratios between the energy densities of the primary 
particles (either CRIs or CREs), magnetic fields, and 
thermal plasma Uth ~ nT. Other possibilities involve a 
primary CRI distribution that has energy density Ui on n 
(such a scaling is less likely for the magnetic fields or the 
rapidly cooling CREs) , or a magnetic field frozen into the 
plasma, B cx n^/'^. The synchrotron emissivity in each 
of the nine model variants corresponding to these three 
primary distributions folded with different magnetization 
levels and scalings, is shown in Table [51 



TABLE 2 

Radio emissivity in model variants with different 
distributions of cres and magnetic fields 



Primary scaling 


CREs 


CRIs 


CRIs 


4- Magnetic scaling 


Uc oc nT 


Ui oc nT 


Ui ocn 


Bcmb > B CC VnT 


jy oc ri^T^ 


jy oc n'^T^ 


> oc n^T 


Bcmb > B oc n2/3 


ju oc n'^/^T 


> oc niO/3T 


> oc 


B > Bcmb 


jv oc nT 


ocri^T 


oc -n? 



Among these model variants, only the two secondary 
CRE models in which the magnetic field is strong (models 
highlighted as boldface in the table) are consistent with 
the scaling Eq. (O — jj. oc n'^T'^ where k = 0.2 ± 0.5 — 
and with Eq. ([IU|. Both models are consistent with our 
data. A slightly better fit to the GH profiles is obtained 
with Ui oc n (i.e. k = 0; see Sj3|), but more data is needed 
in order to establish the thermal dependence of the CRI 
distribution with sufficient statistical significance. 

An independent argument in favor of these two mod- 
els is the environment of MHs. MHs are found in re- 
laxed clusters such as A2029, where only CFs reveal de- 
viations from hydrostatic equilibrium. Such CFs reflect 
subsonic bulk shear flows and strong magnetization, but 
were not associated with particle acceleration. The CF- 
MH connection therefore supports both secondary CREs 
and strong magnetization in MHs. The similarity in the 
values of rj we flnd in GHs and in MHs implies, transi- 
tively, that the same holds for GHs. 

Additional evidence supporting the presence of strong 
magnetic flelds and the absence of primary CREs is found 
by examining the morphological and spectral properties 
of halos, as discussed in Notice, for example, that the 
brightness correlation shown in Fig. ([2]) is strongest in the 
centers of halos, but diminishes at larger radii, where 
turbulent activity associated with mergers is expected 
and where merger shocks and relics are found. In primary 
CRE models, one would not expect the correlation to 
preferentially tighten away from the turbulent regions. 

4.2. Temporal variations 

The preceding discussion is strictly valid only as long 
as the timescale for substantial changes in the magnetic 
held is longer than the cooling time of the CREs, tcooi- 
For CREs that emit synchrotron radiation received with 
characteristic frequency v, 



^cool 



0.13 



(*5) 



-3/2 



^(1 



Gyr (11) 
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where the term in square brackets peaks at unity when 
B = Bcmb/V^, and scales as B~^/^ for B ^ Bcmb- In 
GHs, tcooi is much shorter than the > 1 Gyr timescale 
characteristic of the halo lifetime (K09). 

However, significant local variations in the magnetic 
field can take place on a shorter timescale, of the or- 
der of the sound crossing time of the turbulent eddies, 
tg ^ l/cs- This is shorter than tcooi for a sound veloc- 
ity Cs ~ 10'^ km s~^ and eddy length scales I < 100 kpc. 
(Note that substantial magnetic power is measured on 
coherence scales ~ 10 kpc; see The radio emis- 

sion should therefore be averaged over tcooi and over the 
beam. Nevertheless, as long as many eddies contribute to 
the emission, and the variations in magnetic energy den- 
sity remain of order unity, this correction would be small. 
Notice that fast changes in magnetic configuration over 
a light crossing time, if present, may be observed as tem- 
poral radio variations by next generation telescopes such 
as the Square Kilometre Array (SKA). In deed, the milli- 
arcsecond resolution attainable by SKA (jSchilizzi et al.l 
[2007) corresponds to a light-crossing time of less than a 
year, for nearby halos at redshift z < 0.02. 

In MHs, the magnetic fields are probably associated 
with sloshing activity in the core (see ^6.ip . The charac- 
terist ic timescale for the deca y of core sloshing is > 1 Gyr 
{e.g., lAscasibar fc Markevit ch 2006), much longer than 
tcooi- The timescale for the buildup of sloshing depends 
on its trigger mechanism; in a merger induced scenario 
this is again > 1 Gyr. However, local variations in the 
magnetic field could occur over the radial sound cross- 
ing time, ts ^ r/cs, or on the crossing time of the 
magnetic structures associated for example with CPs. 
These timescales could be shorter than tcooi in the cen- 
tral ^ 50 kpc (note that in these regions B usually sig- 
nificantly exceeds Bcmb', see ^5.11) . As in GHs, averaging 
over tcooi and over the beam could introduce small vari- 
ations in radio brightness, in particular at the edges of 
MHs. 

In both GHs and MHs, Q/n"^ is weakly sensitive to 
variations in CRE injection, through changes in the frac- 
tional energy of the CRIs, Ui/uth- Thus, Q/n'^ should 
be replaced by its value averaged over tcooi and over the 
beam. This correction should be very small, except near 
a cosmic-ray source. 

5. SIGNATURE OF SECONDARY CREs IN STRONG FIELDS 

In the preceding sections we showed that observations 
support a halo model which invokes secondary CREs 
in strong magnetic fields as the origin of radio halos, 
both GHs and MHs. In such a model, the radio emis- 
sivity ji, (X UpuB^/lB"^ + i?cm&) depends weakly on B. 
This model recovers the observed radio- X-ray correla- 
tions, provided that the local energy fraction of CRPs is 
narrowly distributed about a universal, weak function of 
the cluster parameters, because then jj^ ~ ~ jx. 

Such a model entails particular morphological and 
spectral properties of GHs and MHs. Utilizing the model, 
we show how these properties can be used to shed light 
on halo observations, test the model, and gauge its pa- 
rameters. In particular, we demonstrate how the distri- 
butions of CREs and magnetic fields can be measured 
separately, rather than their degenerate product as done 
in most other models. 



5.1. Morphology: B < Bcmb radio suppression 

In our model, both GHs and MHs are regions in which 
strong magnetic fields with B > Bcmb indirectly illumi- 
nate the cluster's CRI population in radio waves. This 
explains the spatial coincidence between MH edges and 
CPs, which are present in more than half of all CCs. Bulk 
shear fiow is believed to magnetize the plasma across and 
beneath the CPs; t here is however no evidence for shear 
above the CPs ( Keshet et al.|[20Tol ). Therefore, observa- 
tions of sharp MH termination coincident with CPs re- 
flect the transition from strong to weak magnetic fields. 
This is probably the reason for the rapid, nearly exponen- 
tial ri{r) cutoff in Perseus above r ~ O.OSrsoo — 50 kpc, 
shown in Pig. [51 Note that CPs are not spherical; the ra- 
dial decay may result from one or several CPs extending 
over various radii, seen projected and radially binned. 

In GHs, and in MHs away from CPs, it is natural to 
assume that the magnetic field decays with r, possibly 
as some fixed fraction of equipartition, B^ cx nT. At 
some distance rf,, B drops below Bcmb, leading to a ^ 
[B / BcmbY' suppression in the radio emissivity. This can 
produce a radial break in the projected radio profile, with 
77 ~ constant at r ^ rf,, and rj (x B^ (x r*^^ at r ^ rj, 
with some < 0. A power-law pressure profile nT cx 
r"^ would imply kb = up', for an isothermal distribution 
= -2. 

Such a decline in rj is found in A2029 above rf, < 
120 kpc, and possibly also in the MHs in A1835 and 
RXJ1347, as seen in Pig. [31 The figure also shows a sim- 
ple model for A2029, where we assume that B^ = knT 
and that rj asymptotes to its average GH value at small 
radii, ri{r — > 0) — > 770 = 10~'*'^, so the proportionality 
constant k is the only free parameter. This fit (dot- 
dashed curve) corresponds to a central magnetic field 
Bq — 2.1 fiG. However, due to the limited range of data 
in this cluster, there is a degeneracy between Bq and tjq 
— higher magnetic fields are possible if the central value 
of rj in A2029 is lower than the GH average, and vice 
versa. The exponential fit of M09 (dashed red curve) 
suggests higher magnetic fields and lower rjo', a fit with 
Bq = 6.3 fj,G is shown (dotted) to illustrate this. The 
corresponding low value of 770 in this cluster could arise 
from its low star formation rate, as discussed in !j6.2l 

Assuming that B^ cx nT for r < ri,, we can place a 
lower limit on the central magnetic field amplitude Bq in 
each observed halo, given a pressure model. If no break 
in 77 is identified out to the halo radius R^, then 



Bo > Bo.n 



{nT)r=Q 

{nT)r=R^ 



1/2 



Bcrnbiz) ,(12) 



where we allowed a fraction fic ~ 0.5 of inverse Comp- 
ton losses at R^. The Bo.min values of the halos in 
our sample are presented in Table [I] assuming that 
Ri, ~ min{rh,3re} and adopting individual isothermal 
/3-models for each cluster, as detailed in the table. In the 
MHs, the central magnetic field could be much higher 
than Bodmin because (i) the significant growth in nT to- 
wards the center (due to the density cusp) is not captured 
by the /?- model; and (ii) the magnetic field immediately 
bene ath the CPs could be very strong, near equiparti- 
tion ([Keshet et al.|[20Tof ). However, in MHs with only a 
declining rj profile observed {e.g., A2029), adopting the 
upper limit on may overestimate Bq. 
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Instead or in addition to a radial break in brightness, a 
halo in which the magnetic field is marginal, B ~ Bcmb, 
could become clumpy or filamentary, appearing bright 
only in B > Bcmb islands. This is probably the reason 
for the unusual clumpy or filamenta ry radio morpholog y 
observed in RXC J2003.5-2323 (Giacintucci et alJl2009[ ). 
A2255, and A2319 (M09). Indeed, a decline in r] is 
directly seen in A2319 around r = 150 kpc (suggest- 
ing low magnetization; see Fig. [2]), and A2255 is the 
only strongly (20%- 40%) polarized GH known todate 
(jGovoni et al.l feOOS): the absence of strong beam depo- 
larization sug gests relative l y wea k magnetic fields (see, 
for example. iMurgia et al.|[2004[ ). Additional evidence 
for the low magnetization in these three halos is their 
relatively steep spectrum, as discussed in 

Such marginally magnetized halos are particularly in- 
teresting because the magnetic field can be determined 
in multiple locations, and because the radio morphology 
directly gauges the magnetization or magnetic decay pro- 
cess. Fast putative variations in the magnetic field may 
be easier to detect in such regions, as they could involve 
temporal changes in the small scale radio morphology. 
For example, in A2319, SKA could resolve such changes 
over a light-crossing time of ~ 2 years. 

While radio emission rapidly declines outside the B > 
Bcmb region, an opposite signature is expected in inverse- 
Compton emission from the same CREs, as they scatter 
CMB photons to high energies. This emission may be 
clumpy or filamen tary in .B ~ Brmh regions. However, 
as pointed out bv lKushnir fc WaxmanI (pOlOf ). such ra- 
diation cannot explain the hard X-ray excess detected 
in a number of clusters (for review, see iRephaeli et al.l 
|2008[ ). because the inverse Compton signal from sec- 
ondary CREs has i^Ii, comparable to that in the radio 
(see Eq. (jl))), ^ 3 orders of magnitude lower than needed 
to account for the detected hard X-ray excess. 

5.2. Using the radio spectrum to disentangle CRIs and 
magnetic fields 

Spectral steepening with increasing r, increasing v, 
or dec reasing T, has been repor t ed in several radio 
halos (iFeretti fc Giova nninil 120081: iFerrari et all 120081: 
iGiovannini et al. 2QQ9. . and references therein). Such 
trends naturally arise in our model, due to the energy- 
dependence of the inelastic cross-section for collisions of 
CRPs with the intracluster gas at the relevant Ep ^ 
20 GeV CRP energies. 

Secondary GREs are mainly produced by charged pion 
production, p + p tt* -|- X, where X is any combina- 
tion of particles, followed by mesonic decays tt* — fi^ 



'^fj,{i^fj.) and leptonic decays /z^ 



Other processes, in particular He-p collisions, have a 
similar cross-section per nucleon and should not mod- 
ify our results by more than 10%. On average, the en- 
ergies of the GRE, the tt*, and the CRP are related 
by Ec ~ /ei Tr-^Tr ^ fe^UlpEp, where /eu ~ 1/4 and 
U\p ~ 1/5 (jGinzburg fc Svrovatsk\l 119611 ). The inclu- 
sive cross section for tt^ production varies significantly 
as a function of Ep around Ep ~ 10 GeV, dropping from 
> 40 mb above 30 GeV to ze ro at the threshold energy 
0.3 GeV (|Blattnig et al.l[2000l) . 

We compute the radio spectrum based on two dif- 
ferent methods. Model A uses the spectral fits for e'^ 



and e production i n ine lastic p-p scattering accord- 
ing to iKamae et all (|2006l vahd for 0.5 GeV < Ep < 
512 TeV), with corrected parameters and cutoffs (T. Ka- 
mae & H. Lee 2010, private communication). Model 
B assumes that Eg = {Ep/20) and utilizes the inclu- 
si ye cross sections for charged pion production according 
to lBlattnig et al.l (jlOOO), which a gree with expe rimental 
data for 0.3 < (Ep/GeY) < 300 (iNorburvl 120091) . Model 
A should be more accurate, because feiTrfnlp is not inde- 
pendent of Ee and the CRP spectrum. 

The results, depicted in Fig. [3] for emitted frequency 
Vc — [1 + z)v — 1.4 GHz and various power-law CRP 
spectra np{Ep) cx E^, show that the radio spectrum 
steepens with increasing cosmic-ray electron/positron 
energy. 



Ec ~ A.2vl^l^{B/5 AiG)"i/2 GeV 



(13) 



Model A features a spectral break at Ee — 2 GeV, at- 
tributed mainly to electron production through diffrac- 
tive processes, in which one or both protons transi- 
tions to an excited state (discrete resonances and con- 
tinuum). Note that even if this channel is blocked, a 
similar spectral break is found at Ec ^ I GeV. A more 
dramatic steepening than in either model is found if 
we assume that Ec = (E-^/A) and use the pion spec- 
tra produced in p-p collisions according to lBlattnig et al.l 
^000); however, these formulae have not been tested be- 
yond Ep ^ 50 GeV. 



-0.7 



X 

o 




EJGeV 

Fig. 3. — The relation between a (dashed curves; left axis), B 
(solid curve, right axis) and Ee (abscissa), for fe = 1.4 GHz. The 
highly magnetized, B^mbi^ = 0) < B < 200 fiG regime is shaded. 
The radio spectral index ai.4 is shown as a function of Ee according 
to model A (dashed) and model B (dash-dotted), for CRP spectral 
indices s = —2, —2.5, and —3 (top to bottom). See ij5.2l for model 
and parameter description. Inset: The CRP energy index s as a 
function of radio spectral ind ex a 1.4 for B = Bemb (for example, 
at the spatial halo break; see i]5.1ll . 

Due to the limited energy range over which the en- 
ergy index Sc = dlog(ne)/(ilog(i?e) of the injected GREs 
can be computed, the radio spectrum in Fig. [3] is calcu- 
lated in the approximation where ea ch CRE emits a sin - 
gle photon, such that a ~ Se/2 f e.g.. 'Kes het et ani2003f) . 
A more accurate computation, convolving the CRE dis- 
tribution with the synchrotron emission function, would 
somewhat smear the spectral features. 

The radio steepening with increasing Ec implies steep- 
ening with increasing v or r, or decreasing T, as ob- 
served. Both models agree that spectral steepening by 
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|Aa^| > 0.2 indicates B > 2{){v^/lA GHz) mag- 
netic fields in the region associated with a flatter radio 
spectrum. Measuring a across the — 2 GeV spectral 
break can be used to unambiguously fix the CRP spec- 
trum s. Similarly, measuring a at the spatial rj break 
where B ~ Bcmb can be used to determine s, as illus- 
trated in the inset of Fig. [31 Conversely, if a distinct 
spectral break exists as predicted by model A, it would 
directly gauge the magnetic field, once the CRP spec- 
trum has been determined. Multi-frequency radio data 
could thus allow a sensitive mapping of the magnetic 
field, both above and below Bcmb, throughout the halo, 
using, for example, deprojected maps at several radio 
frequencies. 

Spectral measurements of radio halos should be inter- 
preted with caution. Contamination by CRE sources 
such as shocks, relics, central AGNs, and radio galax- 
ies, is common. Spatial averaging often blends together 
different sources, depending on flux sensitivity and an- 
gular resolution. Therefore, with present observations, 
the spectrum can be reliably associated with a halo only 
when measured locally in uncontaminated regions; in our 
model these are regions of uniform rj. Also note that ra- 
dio emission above 1 GHz is i ncreasingly s uppressed by 
the Sunyaev-Zel'dovich effect (|Enfilinll200^ . 

For example, ri{r) is nearly constant in A665 within 
the range r < Tmax — 150" examined in Fig. [2j Spectral 
stee pening with increas ing r was identifled in this clus- 
ter (jFeretti et a l."2004b") by comparing radio maps at 0.3 
and 1.4 GHz frequencies. The spectral index steepens 
from a ~ —1.0 in the halo's center to a ~ —1.3 towards 
the South and towards the East, within the flat r] re- 
gion. Similar steepening was found towards the North 
and West beyond Vmax, but the spectrum first fiattens 
to a ~ -0.9. 

Similarly, the regular GH in A2744 exhibits a linear 
radio-X-ray correlati on in brightness out to r ~ 1 Mpc 
(|Govoni et al.|[2001a|) . Along the main NW elongation, 
the 0.3-1.4 GHz spectral slope a{r < 200 kpc) ~ -1.0, 
slightly fiattens to —0.9 around 500 kpc, and steepens 
again to — 1.5 as r — >■ 1 Mpc. An azimuthal average, how- 
ever, includes a NE relic tail and under-threshold regions, 
leading to an unrealistically uniform a{r < Mpc) ~ — 1 
A similar behavior is found in A2219 
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In these examples, the spectral index steepens with 
increasing r from a ~ —1.0 to —(1.2-1.5) in the less per- 
turbed/contaminated direction, sometimes after a mild 
flattening to a ~ —0.9. Such steepe ning was also 
reported in A2163 (libretti et all [2004bD and in Coma 
(|Giovannini et al.lll993[ ). 

Comparable steepening was found as a function of 
frequency in the GH in A754. Its spectrum measured 
between 74 and 330 MHz, a'^ Q-^ ~ —1-1, steepens to 
aj;! ^ —1.5 (jBacchi et al.ll2003[) . More substantial steep- 
ening has b een reported in other clusters, such as A231 9 
(jFeretti et al. 1997) and A3562 ( Giacintu cci et al.|[2005l ). 
However, contamination by extended radio galaxies was 
reported in these halos. 

Significant steepening with increasing r or v as in the 
examples above is more consistent with the ^ 2 GeV 
spectral break of model A than with model B. In model 
A, strong magnetic fields with B > 10(i/e/700 MHz) ^G 



are present in regions where is fiat. Moreover, an un- 
contaminated spectral break measured at some emitted 
frequency Vf, would imply that the local projected mag- 
netic field is i? ~ 10(i^b/700 MHz) /iG. Observations at 
several frequencies can thus be used to map the local 
projected magnetic field. As a preliminary example, in- 
terpreting the radio spectrum of A754 as a spectral break 
somewhere between i/g — 200 and 900 MHz near the cen- 
ter of the halo, would imply 3 < Bq/ < 14. 

Figure [3] indicates that radio steepening from a = —1.0 
to —1.3 corresponds to a CRP spectral index s ~ —2.7. 
Steepening beyond —1.5, if uncontaminated, would im- 
ply a steep CRP spectrum with s < —3. The radio spec- 
trum at high, Ee » 10 GeV CRE energies approaches 
a — s/2. However, inward of the spectral break, at 
Ee < 2 GeV, a ~ — 1 depends only weakly on s. This 
could explain the universally observed a ~ —1 in the 
centers of halos. 

Figure m shows the average i/ < 1.4 GHz spectral in- 
dices (ai.4) reported for GHs in our sample, as a function 
of their minimal central magnetic fields -Bo.min calculated 
in i 35.ll The results are also summarized in Table [T] A 
positive correlation is seen, in the sense that halos with 
stronger central magnetization tend to have a fiatter ra- 
dio spectrum, as expected. Note that this correlation is 
stronger than, and not simply due to, a correlation we 
find between (a) and the halo size. 

For comparison, we also plot (dashed lines in Fig. HJ 
the average radio spectrum (ai.4) as a function of Bq, 
computed for a typical isothermal /3-model with /3 = 0.7, 
assuming B^ oc nT, spectral model A, and radio emis- 
sion extending to Ri, ~ Sr^- Although this result cannot 
be quantitatively compared to the GH data (which are 
only lower limits on Bq), the trends qualitatively agree. 
This provides an independent indication that the CRP 
spectrum is steep, with s < —2.7. 
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Fig. 4. — The relation between the spatially averaged, 1/ < 
1.4 GHz spectral index (ai.4), and the central magnetic field Bq. 
The measured spectral index in each GH is plotted as a f unct ion of 
the minimal central magnetic field -Bo.min calculated in JSTT] (right 
arrows; data with no reported error bars shown as diamonds; data 
with error bars larger than the y-axis shown as a circle; see Table 
[TJ . Also shown (dashed curves) is the average spectrum calcu lated 
in a simple magnetic equipartition, isothermal /3-model (see i|5.2ll . 
as a function of the central magnetic field, for CRP energy indices 
s = —2, -2.5, and -3 (top to bottom). 



The radio spectrum is relatively steep in the five GHs 
and MHs for which we presented evidence for low mag- 
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netization in ^5.11 in agreement with the (a) — Bq cor- 
relation. The halos which are not shown in Fig. [4l 
include A2029 ((ao.2) ^ -1. 35, [leej^J iegmah logs'), 
A2255 {{ali) = -2.06±0.57. IKemDner fc Sarazin 2001, 
where we inf er {a })'f,) ~ —1.5 by comparing MOO with 
^ " ^' -1.2.! Gitti et al.. 

- -1.32 ± 0.06, 



Harris et all (|1980D V Perseus ((aj; 



200l. and RXC J2003 .5-2323 ((aj-^ 
Giacintucci et al.l l2009V 

Finally, a cautionary comment is in place regarding 
the extrapolation of the spectrum to low frequencies 
that correspond to i?e ^ 1 GeV ORE energies. Al- 
though Ee = 1 GeV CREs are produced on average by 
Ep ~ 20 GeV CRPs (for s = -2.5), the contribution of 
low energy CRPs with Ep <3 GeV is not negligible. We 
expect a CRP power law energy spectrum np{Ep) (x E^^ 
to be modified around the proton rest mass. Conse- 
quently, our s = constant approximation becomes in- 
creasingly unlikely at low, E^ ^ few GeV CRE energies. 

6. DISCUSSION: ORIGIN OF CRIs AND MAGNETIC FIELDS 

In the previous sections we studied the luminosity (in 
fj^l) and the surface brightness (in @ properties of GHs 
and MHs, showed that they are best explained by a 
model invoking secondary CREs and strong magnetic 
fields as the source of radio halos (in 2]), and discussed 
some implications of such a model (in fj5l- Here we con- 
sider the origins of the magnetic fields ft j6.ip and of the 
primary CRIs ft ^6.2p . We subsequently summarize our 
work in t ^6.3l 

6.1. Magnetization: mergers (GHs) and sloshing (MHs) 

In the model derived above, the primary CRI popula- 
tion is long-lived and has very similar properties in dif- 
ferent clusters and in different parts of a cluster. Hence, 
the defining property of halos is sufficiently high magne- 
tization. Measuring strong, B > Bcmb magnetic fields 
without radio detection at the level given by Eq. Q 
would rule out the model, unless the CRI energy frac- 
tion is exceptionally low (e.g., due to a low level of star 
formation; sec HQ.2\i . The model can be tested by check- 
ing if the different magnetic field estimates it produces 
in a given halo are self consistent, and agree with inde- 
pendent measurements. 

In the model, GH clusters are associated with strong, 
B ^ Bcmb ^ 3 /iC magnetic fields, whereas weaker, 
B < Bcmb fields are present in clusters devoid of a 
halo. Magnetic fields of the order of Bcmt, extending 
over Mpc scales, are consistent with Faraday rotation 
measures (RM) in non-cool clusters, considering the un- 
certainties involved. Such RM studies suggest magnetic 
fields ranging between a few /iC to 10 /iC, extending out 
to r ~ 500 kpc, in a sample of non-cool c ore clusters, on 
10-20 kpc coherence scales (jClarkell2004[) . Weaker mag- 
netic fields, typically ^ 1 /xG to a few //G, were found 
using other methods, such as simulating the polariza- 
tion of extended sources (Murgia et al. 2004) or using 
a Bayesian maximum likelihood analysis of the Faraday 
RM ( Vog t k Enfilin 2005). 

RM studies could potentially test the association be- 
tween strong magnetization and the presence of halos, 
as implied by the model, although this is complicated 
by several factors. First, RM analyses involve substan- 
tial systematic and statistical errors associated with pro- 



jection effects, the magnetic power spectrum, the range 
of coherence scales, and non-Gaussianity. Second, clus- 
ter fields are often reconstructed using peripheral RM 
sources, assuming some magnetic field scaling B on n'^ 
that extends both inside and outside the halo; however 
this assumption is uncertain and in some cases incon- 
sistent with our model, e.g., in asymmetric or clumpy 
halos. Finally, the RM gauges the magnetic field ampli- 
tude integrated along the line of sight, whereas the radio 
emission scales differently with the magnetic field, e.g., 
as B^ in B < Bcmb regions. 

It is interesting to compare the magnetic fields es- 
timated in different clusters. Care must be taken to 
compare the same measures estimated under the same 
assumptions. For example, assuming a Gaussian, Kol- 
mogorov magnetic power spectrum with fj, = 1/2, we 
find ( B(r < Src)) = 1.8 uG and B{3rc) = 1.3 fiG in 
Coma (jBonafede et al.ir2010( ). which harbors a GH, while 
weaker fields, {B{r < 3rc)) = 1.2 /xG and Bj Src) = 
0.8 AiG are found in A2382 (jGuidetti et al.ll2008[) . which 
has no halo. Similarly, {B{r < 3 rr)) = 1-3 mG an d 
B{3rc) = 0.8 A*G are found in A119 C Murgia et a"l1 l2004V 
which does not harbor a GH, but this result is obtained 
for /i = 0.9 so comparison to the other two clusters may 
be misleading. Note that it is not clear, considering the 
inherent uncertainties, if the differences between these B 
estimates are significant. 

The association between GHs and merger events is 
fairly well established (see iFeretti Sz Giovanninil I2008L 
and references therein). With present magnetic field es- 
timates, it is quite plausible that a major merger event 
could amplify the intracluster magnetic field by a factor 
of a few, sufficient to exceed Bcmb and produce a GH. 
The model thus resolves the notorious discrepancy be- 
tween magnetic field estimates based on Faraday RM and 
on previous GH analyses, which unnecessarily assumed 
equipartition between magnetic fields and CREs. 

The very strong magnetic fields we infer in the centers 
of MHs are consistent with the RM observed in CCs, 
which typically imply B ~ 10 — 40 u G in the core, on 
coher ence scales of a few up to 10 kpc (jCarilli fc Tavloil 
|2002| ). It is interesting to point out that the cooling 
flow p ower is correlated both with the RM (iTavlor et aLl 
120021) and with the MH radio power (jGitti et al.ll2004D . 
while no strong corr elation between MH and AGN power 
was identified (e.g.. iGovoni et al7 '2009\ 

The association between MH edges and CFs strongly 
suggests that sloshing motions play a major role in mag- 
netizing the core. Notice that ~ 10 kpc scales are in- 
deed characteristic of the shear magnetic am plification 
anticipated beneath CFs (iKeshet et al.l^ 20Tol ). Sloshing 
could suppress cooling in the core , for example by mixing 
the c old gas with a heat inflow ( Markevitch fc VikhlininI 
[2OOI . In such a scenario, a stronger cooling flow may 
correspond to a larger magnetized region, leading to cor- 
relations between the cooling flow power and both RM 
and MH power, as observed, but not necessarily to an 
MH-AGN correlation. 

The distribution of MHs (without the AGN compo- 
nent) among CCs should reflect the statistics of core 
magnetization, and therefore of sloshing. We expect a 
correlation between the presence of MHs and of CFs in a 
cluster, and a correlation between the MH size or power 
and the shear flow strength, manifest for example in the 
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number and size of CFs and in the (measurable, see 
iKeshet et al.|[20lol ) shear across them. Assuming that 
some level of magnetization by core sloshing is always 
present in CCs, as suggested by the ubiquity of observed 
CFs, the steep magnetic profile associated with the den- 
sity cusp would imply that every CC has some MH, even 
if small. The MH distribution would then be continuous, 
and not bimodal as in GHs. 

Shear amplification of magnetic fields parallel to the 
CP plane is expected mostly below the CF and in a thin 
boundary layer around it, in which the field can reach 
high, near equipartition levels (jKeshet et al.ll201Clf ). The 
morphological association between MH edges and CFs, 
discovered by iMazzotta fc Giacintuccil (|2008^ , is there- 
fore expected to be ubiquitous. It is best seen where 
CFs are observed edge on; elsewhere it may be observed 
as morphological correlations between radio maps and 
spatial X-ray gradients, which trace projected CFs. 

The magnetic field (the polarization) is expected to be 
parallel (perpendicular) to the CF, i.e. approximately 
tangential (radial). Polarization would be preferentially 
observable where beam deprojection is minimal, i.e. at 
large radii where the magnetic field weakens. Interest- 
ingly, nearly radial, 10% — 20% polarization was de- 
tected in the MH in A2390, growing stronger with in- 
creasing radius (we refer to the spherical, r '-^ 150 kpc 
co mponent around th e cD galaxy in this irregular MH; 
see lBacchi et al.ll2003| ). We predict that near CF edges, 
where the magnetic field is particularly strong, the radio 
spectrum would be relatively fiat (see §5.2|) . 

6.2. CRI Origin: diffusion and supernovae sources 
The spectral steepening of the radio signal with in- 
creasing CRE energy Ee oc {v/By^^, provides a novel 
method for measuring the primary CRP spectrum. The 
radio steepening observed in some halos, roughly from 
a ~ —1.0 to a ~ —1.3, corresponds to a CRP spectral 
index s < -2.7 at Ep - 20 GeV energies. The CRP 
energy fraction then becomes (c/. Eg. [TU| 

f (-^ r \ — '"p(> ^ in-3.6±0.2 { \ (-,A\ 

Note that with the uncertain and possibly contaminated 
radio spectra presently available, a steeper CRP spec- 
trum with s ~ —3 is possible. The model would be 
challenged if the uncontaminated spectrum of a substan- 
tial halo population turns out to be much steeper than 
a = —1.5, unless the corresponding steep CRP spectrum 
can be explained. 

The CRP distribution in Eq. resembles (but 

has an energy fraction a few 100 times smaller than) 
the CRP distribution found in the solar vicinity above 
1 GeV/nucleon. This distribution could originate from 
sources that inject roughly equal energy per decade of 
CRP energy (sq — —2.2), such as supernovae (SNe), 
if energy-dependent diffusion is significant in the in- 
ner halo regions. For example, a simple estimate of 
CRI scattering off magnetic irregularities with a Kol- 
mogorov power spectrum yields a diffusion coefficient 
D ~ 1030(£;p/GeV)i/3(B///G)-i/3 cm^ s"! (|V61k et al.l 
I1996D . This imphes CRI diffusion over ~ 0.5 Mpc during 
a Hubble time, and a steepening by As = —1/2. More 
substantial steepening is possible if the diffusion function 



has a stronger energy dependence. For example, the dif- 
fusion function is often assumed to scale as c>c E^^^ , 
which could lead to a As = — 3/4 steepening in the CRI 
spectrum. 

The CRI outp ut of SNe can be crudely estimated 
(IVolk et al.lll996l) if we assume that a fraction fu of the 
cluster's Z = 0.3.^0.3 solar metaUicity is seeded by Type 
II SNe, which on average produce 0.1MQMpe,o.i of iron 
and deposit a fraction ^ = 0.3^o.3 of the lO^^i^si erg ex- 
plosion energy in £'p > 10 GeV CRIs, 

~ 0.03/iiZo.3£;5lMf e,o.iCo.3 ■ (15) 

This can reproduce Eq. if over the cluster's lifetime, 
the CRIs diffuse to distances a few times larger than the 
radius Ri, of the radio halo. Note that if CRI diffusion 
is entirely absent, the CRIs accelerated in SNe would be 
confined to the cluster, and adiabatic losses could only 
lower their energy density to the level of Eq. ([T^ . How- 
ever, they would then retain their flat, s ~ —2.2 spec- 
trum. 

An SNe origin of CRIs can be tested by examining the 
correlations between rj and (intensive) tracers of SNe ac- 
tivity among different halos. One possible tracer is the lo- 
cal metaUicity measured at r = O.lrsoo, tabulated in Ta- 
ble[TJ We chose to use Z(0.1r5oo) because (i) it was mea- 
sured for all the M09 halos with XMM-Newton profiles in 
ISnowden et~aLl (|2008 f): /mJ the spatially averaged Z is not 
meaningful when temperature gradients are large; and 
(Hi) O.lrsoo lies well within the Z{r > 0.02r5oo) c>c r~°-^ 
decline t ypically found in both cool and non-cool core 
clusters ([Sanderson et al.l I2009D . While our sample is 
statistically small, Perseus, which shows a significantly 
higher 77 than in all the other hal os in our sample, als o 
shows a slightly higher ^(O.lrsoo) (|Snowden et al.l l2008V 
However, due to the large uncertainty in abundance mea- 
surements, the elevated metaUicity in Perseus is not sig- 
nificant (< la) with respect to some halos. Moreover, 
at smaller radii r < 0.025r5oo, the metaUicity in A2029 
appears to be higher than in all other halos, and is sig- 
nificantly (> 3cr) higher than in Perseus (Snowdcn ct ah] 
1 20081) . 

Although better metaUicity statistics may identify a 
more significant correlation between rj and Z, metaUic- 
ity is probably not the most useful tracer of the SNe 
contribution to CRIs in the halo. MetaUicity provides a 
cumulative measure of SNe activity, tracing the metals 
released from all past SNe in the cluster. In a model 
where a significant fraction of the CRIs have already dif- 
fused away from the cluster's center, metaUicity would 
not linearly trace the population of CRIs residing within 
the halo, especially in the more compact MHs. It is more 
appropriate to use an intensive tracer of recent SNe ac- 
tivity, such as the star formation rate (SFR) normalized 
by the cluster's gas mass Mg, or the fraction of star form- 
ing galaxies. A correlation between an SNe measure and 
a CRI tracer, such as t] or the deviation from the lu- 
minosity correlation i^Piy/L^J^ 124]' ^® more useful 
than the metaUicity in establishing or ruling out an SNe 
origin of halo CRIs. 

As seen in Fig. [51 the halos with the highest 77 in our 
sample are the MH in Perseus and the GH in A665, while 
the halos with the lowest rj are the MH in A2029 and 
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the GH in A773. Interestingly, the hterature shows ev- 
idence for exceptionally high specific star formation in 
both Perseus and A665, and for a low specific SFR in 
A2029. (We found no relevant data for A773.) 

In Perseus, which has the smallest Mg (by at least 
a factor of 4, see .Fukazawa et al.. .2004i) and one of the 
most power f ul cooling flows w ithin our sample {e.g., 
iWhitd [2OOOI: lAllen et al.|[200^ . there i s optical-to-UV 
evidence for a relatively h igh SFR {e.g., iBregman et alj 
l2006HRaffertv et al.l[2008l ). In particular, the cD galaxy 
NGC1275 in Perse us has a high SFR of ~ 30Mq yr"! 
([Dixon et al.l 119961 ) — the highest in our sample when 
normalized by Mg. Notice that the central galaxy in 
A1835 has a higher SFR of ~ lOOMfr, yr"^ — the highes t 
SFR known in such objects (jPeterson fc FabianI 120061 ) . 
However, Mg is ~ 10 times larger in A1835 with respect 
to Perseus (jPukazawa et al.l I2004D . Note that regions 
containing a high density of cosmic-rays are directly ob- 
served in Perseus in the form of X-ray cavities, reaching 
distances r > 100 kpc (M. Markevitch, private commu- 
nications). A high specific SFR is also inferred in A665, 
which was found to be the cluster with the highest dis- 
persion in color magnitude relation (a known tracer of 
star formation) in a sam ple of 57 X-ray bright clusters 
(|L6pez-Cruz et al.l [20041 ). In contrast, A2029 has a low 
specific SFR, as it was shown to have a SFR ^ 70 times 
lower than in A1835 fflick s k Mushotzkv '2005) while 
its ga s mass is only ~ 1.8 times smaller ( Fukazawa et al.l 
I2OOI . 

While these trends support an SNe origin of CRIs, 
more work is needed in order to quantify their signif- 
icance and compile a comparative statistical analysis. 
Note that the combination of a strong correlation of r/ 
with the specific SFR and a poor correlation with mctal- 
licity, if established, would directly imply that CRI dif- 
fusion is significant. Indeed, it is difficult to explain the 
steep, s < —2.7 spectrum without invoking CRI diffu- 
sion. 

The above estimates of diffusive steepening assume 
that most of the CRIs produced by the sources presently 
dominating the halo have already escaped beyond it. 
This is consistent with the typical SFR peak at z ~ 1, 
and with the above estimates of the halo CRI abun- 
dance and the total CRI output of SNe (c/. Eqs. IT4l 
and I15|). However, such substantial diffusion would in- 
troduce some scatter in the radio-X-ray correlation, de- 
pending on the CRI production history of each cluster. 
Quantitative estimates of the SNe history of GH clusters, 
needed to compute this scatter, are beyond the scope of 
this work. 

6.3. Summary and Conclusions 

We have shown that the radio- X-ray correlation in GH 
luminosity (Eq. [T|) can be generalized (Eq. [5]) to hold 
for both GHs and MHs (Fig. [1]), by correcting for the 
halo size. A universal, linear relation between the ra- 
dio and X-ray surface brightness, 77 = 10~'* '**°'^, was 
presented (Eq. [4] and Fig. [2]). This, combined with the 
radial 77 and T profiles, implies a universal radio emis- 
sivity vj^ = 10-3i-4±0-2„2^2(r/To)°-2±0-5 ergs-i cm'^ 
(Eq. [7]) near the center of halos. We argued that these 
results and their applicability to GHs and MHs alike, 
strongly support one model for all halos, involving sec- 



ondary CREs (injected according to Eq. [TU|) and strong 
magnetic fields with B > Bcmb, while disfavoring other 
models (Table [l]). 

This model makes useful predictions without requir- 
ing additional assumptions or fine tuning. Radio emis- 
sion rapidly fades in regions where B drops below Bcmb, 
producing a distinct radial break {e.g., in A2029 and in 
Perseus; Fig [2]) or a clumpy/filamentary radio morphol- 
ogy {e.g., in RXC J2003. 5-2323, A2255, and A2319) that 
can be used to map B ~ Bcmb contours. Marginally 
magnetized regions with B < Bcmb oc (1 -I- z)^ are char- 
acterized by relatively high polarization and a steeper 
radio spectrum; their morphology traces the magnetic 
evolution and can potentially reveal a temporal signal. 
We expect a higher incidence rate of such transition re- 
gions at higher redshift, while no halos should exist at 
very high redshift. 

Another direct consequence of the model is radio spec- 
tral steepening with increasing E"^ oc ly/B (Eq. IT^ . 
i.e. with increasing r or u ot decreasing T, as in- 
deed is observed. Such steepening, and in particular a 
B ~ 10(:^e/700 MHz) ^G spectral break (Fig. E]), gauges 
the magnetic field, roughly producing an additional B 
contour for each radio map frequency. The spectral 
break could be used to accurately map B throughout the 
halo, using future radio telescopes such as the Murchison 
Widefield Array (MWA0), the LOw Frequency ARray 
(LOFAB0), and SKA. 

A pressure model can be used to extrapolate B 
throughout the cluster. This indicates central magnetic 
fields Bq (Eq. rT2|) that exceed 10 fj,G in most halos (see 
Table [T] for lower limits -Bo.mm)- A correlation between 
the average radio spectral index (a) and Bq, implied by 
the model, was identified in GH data (Fig. S]). 

In our model, any source of strong {B > Bcmb), 
persistent magnetic fields in the intracluster medium 
would have similar properties to radio halos, as long as 
it does not significantly inject additional cosmic rays. 
This may include some extended radio galaxies, which 
were recently found to exhibit p roperties similar to halos 
(|Rudnick fc LemmermanI 120091 ). Conversely, the univer- 
sal value of r] we predict for any highly magnetized, un- 
contaminated region in the intracluster medium provides 
a powerful test of the model. 

The spectral steepening of the radio signal, the univer- 
sality of a ~ —1 in the centers of halos, and the corre- 
lation between (a) and Bo^min (Fig. H]), indicate a steep 
CRI spectrum, s < —2.7, and thus favor significant CRI 
diffusion. In a diffusion model, the most plausible source 
of the CRIs (Eq. [HI) is SNe {e.g., Eq. [TS]). We show 
(in i i6.2|) preliminary evidence for a correlation between 
rj and the SFR normalized by the gas mass Mg, sup- 
porting an SNe CRI origin. None of these properties 
is expected in an alternative model (K09), in which the 
secondary CREs arise from ^ 1 GeV CRPs, which are 
accelerated in the cluster's virial shock and advected in- 
ward with the flow, thus being compressed to ^ 10 GeV 
energies. Note that the data slightly favors a ji, cx ri^T" 
scaling within each cluster (see which is natural if 
CRIs originate in SNe, rather than the ji, cx n^T^ behav- 
ior anticipated if they are accelerated in the virial shock. 

^ http://www.mwatelescope.org 
^ http://www.lofar.org 
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Also note that adiabatic compression of CRIs produced 
at the virial shock and advected with the gas would lead 
to a radially increasing, rj (x n~^^^ profile (due the soft 
equation of state of relativistic particles), which is not 
observed (see Fig.[2|). 

We stress that although our model and the model of 
K09 disagree regarding the origin of CRIs, the CRI dis- 
tribution, the spectral properties of halos, and the role of 
diffusion, we reach the same conclusions regarding the ra- 
dio mechanism: emission from secondary CREs in strong 
magnetic fields. This conclusion is based on (i) the tight 
radio-X-ray correlation in total GH luminosity and the 
GH bimodality (K09); (ii) the tight radio- X-ray corre- 
lation in both coincident luminosity and surface bright- 
ness, in both GHs and MHs, despite their different phys- 
ical properties; (Hi) the strong magnetic fields inferred 
from Faraday RMs in MHs and possibly (see H6.ip also in 
GHs; (iv) the tightening of the brightness correlation at 
small radii, away from merger shocks, radio relics, and 
their associated turbulence; (v) the ji, oc ri^T'^ scaling 
of the radio emissivity within each halo, where k < 1; 
(vi) the coincidence between MH edges and CPs, mani- 
fest as a sharp radial cutoff in rj {e.g., in Perseus); (vii) 
a power-law radial break where 77 (r <ti ri,) ~ constant 
and ri{r ^ rf,) cx B^, possibly seen in the MHs in A2029, 
A1835 and RXJ1347, and in the GH in A2319; (viii) the 
clumpy/filamentary morphology of some halos, where in- 
dependent evidence for low magnetization is present; (ix) 
the spectral steepening and the correlation between (a) 
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and Bq (this suggests strong magnetization provided that 
the CRI spectrum is steep, s < —2.7). 

Each aspect of our model can be tested in the near fu- 
ture. The association between the presence of halos and 
strong, B > Bcmb magnetic fields can be directly tested 
by comparing the magnetization levels independently es- 
timated in halo and in non-halo clusters, as illustrated in 
The secondary origin of the halos can be tested if 
the CRIs are detected directly through their 7r° produc- 
tion; for example, such a detection of a CRI component 
substantially stronger than in Eq. ([T4l) would rule out 
our model. The SNe origin of the CRIs can be tested by 
carefully examining the correlation between a CRI mea- 
sure such as r], and an intensive SNe tracer such as the 
specific SPR, in a sample of halo clusters. 
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